Evolution of 7r° suppression in Au+Au collisions from ^/s^ = 39 to 200 GeV 
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Neutral-pion, 7r°, spectra were measured at midrapidity {\y\ < 0.35) in Au+Au collisions at y/s^-^ 
= 39 and 62.4 GeV and compared to earlier measurements at 200 GeV in the 1 < pr < lOGeV/c 
transverse- momentum (pt) range. The high-pr tail is well described by a power law in all cases and 
the powers decrease significantly with decreasing center-of-mass energy. The change of powers is very 
similar to that observed in the corresponding p+p-coUision spectra. The nuclear-modification factors 
(Raa) show significant suppression and a distinct energy dependence at moderate pr in central 
collisions. At high pr, Raa is similar for 62.4 and 200 GeV at all centralities. Perturbative-quantum- 
chromodynamics calculations that describe Raa well at 200 GeV, fail to describe the 39 GeV data, 
raising the possibility that the relative importance of initial-state effects and soft processes increases 
at lower energies. A conclusion that the region where hard processes are dominant is reached only 
at higher pr, is also supported by the xt dependence of the a;T-scaling power-law exponent. 

PACS numbers: 25.75.Dw 



Large transverse-momentum (px) particles produced 
in high-energy nucleus-nucleus (AB) collisions play a cru- 
cial role in studying the properties of the medium cre- 
ated in relativistic heavy-ion collisions. Most hadrons at 
sufhciently high px are fragmentation products of hard- 
scattered partons and their production rate in vacuum, 
as measured in p+p collisions, is well described by per- 
turbative quantum chromodynamics (pQCD) [ij. In the 
absence of any nuclear effects the production rate in rel- 
ativistic heavy-ion collisions in the pQCD regime, i.e. at 
sufficiently high px , would scale with the increased prob- 
ability that a hard scattering occurs, due to the large 
number of nucleons. This probability is characterized by 
the nuclear thickness function Tab 0- However, such 
scaling has been violated to various degrees depending 
both on collision energy, ^s„„ , and hadron px- At lower 
collision energies, the hadron yield is enhanced above the 
expected scaling. This was first observed in p+A and 
this enhancement is generallv attributed to multiple soft 
scattering ("Cronin effect" Q), and is presumed to occur 
in ion-ion collisions as well. Initial parton distribution 
functions in nuclei (nPDF) are different from those in 
protons 01 • 

Finally, if a dense, colored medium is formed in the AB 
collision, the hard-scattered parton may traverse some 
of it, losing energy in the process. Therefore, the ob- 
served yield at a given (high) px will be lower than 
that expected from Tab scaling, exhibiting "suppression" 
or "jet quenching," described in terms of the nuclear- 
modification factor, i?AA (see Eq. (1)). Alternatively, 
other studies divide the yields for heavy-ion collisions at 
one energy with those for the same colliding species at a 
lower energy Au-|-Au , rather than scaled P+P reference 
data, to study energy and centrality scaling [5|. 

One of the first discoveries at the Relativistic Heavy- 
Ion Collider (RHIC) was a very large hadron suppression 
at high Px (above «3 GeV/c) in = 130 and 200 GeV 
Au-I- Au collisions ^-^] . This suppression was attributed 
to the dominance of parton energy loss in the medium, i.e. 
to final state effects. To test this hypothesis, the same 



measurements were performed in c?-|-Au collisions [lo| . 
where the formation of the hot, dense partonic medium 
is not expected, and initial-state effects (if any) prevail. 
No suppression in d-l-Au data was observed leaving little 
(if any) room for the initial-state effects as the origin of 
the large jet quenching observed in Au -I- Au. Studies 
with the lighter Cu-|-Cu system at three energies (y^s^ 
= 22.4, 62.4 and 200 GeV ^) have revealed that at 



22.4 GeV mechanisms that enhance i?AA (> 1) 
dominate at all centralities. Note, however, that this 
data set had very limited px range {px <4GeV/c). At 
62.4 GeV, jet quenching overwhelmes any enhancement 
and leads to a suppression (Raa < 1) in more central 
collisions. 

The low-energy scan at RHIC provides an opportunity 
to study the transition from enhancement {Raa > 1) to 
suppression {Raa < 1) and the evolution of i?AA with 
collision energy, centrality and px- The results put con- 
straints on energy-loss models (see I4I and references 
therein) . 

Here, we present new measurements by the PHENIX 
experiment at RHIC of tt" invariant yields and i?AA in 
Au + Au collisions at ^/s-^^ = 39 and 62.4 GeV. The data 
were taken during the 2010 run and the px limits (statis- 
tics) were 8 GeV/c and 10 GeV/c, respectively. Reference 
p-|-p-collision data for y/s^ = 62.4 GeV were taken in the 
same experiment in the 2006 run [l|, while for ^/s^ = 
39 GeV, data measured in the FERMILAB experiment 
E706 were used ^1>^. 

Neutral pions were measured on a statistical basis via 
their tt^ — ^ 77 decay branch with the electromagnetic 
calorimeter (EMCal) The EMCal comprises two 

calorimeter types: 6 sectors of lead scintillator sampling 
calorimeter (PbSc) and 2 sectors of lead glass Cerenkov 
calorimeter (PbGl) . Each sector is located w 5 m from 
the beamline and subtends |?7| < 0.35 in pseudorapidity 
and Acf) = 22.5° in azimuth. This Letter presents results 
obtained with the PbSc sectors only. The segmentation 
of the PbSc (A77 xA<j)^ 0.01 x 0.01) ensures that the two 
photons from the vr" — > 77 decays are very well resolved 
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up to pt < 12GeV/c, i.e. across the entire px range of 
this measurement. 

The results are based on data sets of 3.5 • 10® and 7.0 • 
10® minimum bias Au + Au events at 39 and 62.4 GeV, 
respectively. The minimum bias (MB) trigger for both 
y/s^ = 39 and 62.4 GeV was provided by the Beam- 
Beam-Counters (BBC) [l3|, located close to the beam 
axis in both directions and covering 3.0 < \r]\ < 3.9. In 
order to reduce background at least two hits were re- 
quired in both BBC's. This condition selects ^ 86% of 
the total inelastic cross section. The centrality selection 
in Au + Au collisions at both energies was based on the 
charged signal sum of the BBC's, which is proportional 
to the charged particle multiplicity. For each centrality 
the average number of binary collisions ((iVcoii)) and the 
number of participants ((iVpart)) were calculated using a 
Glauber model 2] based Monte Carlo code. 



10 




AuAuVs = 39 GeV,7t" 

G Min.bias 

• 0-10 %x 10'^ 

■ 10-20 %x 10"^ 

A 20-40 %x 10"^ 

T 40-60 %x 10^ 

O 60-86 %x 10"^ 




AuAuVs = 62.4 

□ Min.bias 

• 0-10 %x II 

■ 10-20 %x 

A 20-40 % X 

T 40-60 % X 

yn O 60-86 "Ax 



~ (b) 



O 



T 



P, [GeV/c] 



FIG. 1: (Color online) Invariant yields of in Au -|- Au at 
^/s^ = 39 GeV (a) and 62.4 GeV (b) in all centralities and 
minimum bias. Only statistical uncertainties are shown. 



ing the nuclear modification factor 



TABLE I: Sources of systematic uncertainties and their rela- 
tive effect (in %) on the invariant yields 



-Raa(pt) 



Pt 


2GeV/c 


5GeV/c 


Type 


energy 




39 GeV (62.4 GeV) 




Yield extraction 


3% 


(3%) 


3% 


(3%) 


A 


FID efficiency 


4.5% 


(4.5%) 


4.5% 


(4.5%) 


B 


Energy scale 


10.5% 


(8.0%) 


14.5% 


(10.0%) 


B 


Acceptance 


2% 


(2%) 


2% 


(2%) 


B 


Conversion 


4% 


(4%) 


4% 


(4%) 


B 


Off vertex 


1.5% 


(1.5%) 


1.5% 


(1.5%) 


C 


Total for 7r° yields 


12.7% 


(10.7%) 


16.2% 


(12.3%) 





The PHENIX analysis of neutral pious is described in 
detail elsewhere [9f . TableUlists the sources of systematic 
uncertainties on the extracted- tt*' invariant yields in this 
analysis. They can be divided into three different cat- 
egories: (1) Type-A, pr-uncorrelated; (2) Type-B, px- 
correlated, where the correlation may be an arbitrary 
smooth function; (3) Type-C, pr-correlated, where all 
points move by the same fraction up or down. The main 
sources of systematic uncertainties in the tt*^ measure- 
ment are the energy scale, yield extraction and particle- 
identification (FID) efficiency correction. 

Figure [1] shows the invariant yields of the tt^s for all 
centralities and also in minimum bias collisions. From 
fitting the i/Vn = 39 and 62.4 GeV minimum bias spec- 
tra with a power law function (cx pj.) for pT > 4 GeV/c, 
we obtained powers 71,39 = —13.04 ± 0.08 and ne2A = 
— 10.60 ± 0.03, respectively, significantly steeper than 
at =200 GeV, where nzoo = -8.06 ± 0.01 for 

MB collisions |9j . The slopes of the corresponding p+p- 
coUision spectra are somewhat different, but compara- 
ble, n?g = -13.59 ± 0.21, ngf 4 = -9.82 ± 0.18 and 



'200 



'39 — -L«J-'J» -L u.^±, /tg2.4 

-8.22 ± 0.09, respectively. 



(l/VrA)d^^lI/dpTdy 
(Tab) X d^afJdpTdy ' 



(1) 



Nuclear effects on the tt^ production are quantified us- 



where ct^^ is the production cross section of tt^ in p+p 
collisions, and (Tab) — (Vcoii) /Cpp^' is the nuclear thick- 
ness function averaged over the range of impact param- 
eters contributing to the given centrality class according 
to the Glauber model. Thus -Raa compares the yield 
observed in A + A collisions to the yield expected from 
the superposition of Vcoii independent p+p interactions. 
In the absence of nuclear effects, -Raa should be equal 
to unity. However, -Raa ~ 1 does not necessarily imply 
the absence of suppression, it may also indicate a balance 
between enhancing and depleting mechanisms. 

In order to calculate -Raa, a reference pt distribution 
in p+p collisions is needed. Preferably this is measured 
with the same detector, in which case many systematic 
uncertainties cancel in the ratio. The PHENIX exper- 
iment has measured the tt*^ cross section in p+p colli- 
sions at ^/s^ — 62.4 GeV 1] but only up to pr = 
7GeV/c while the current Au-|-Au measurement reaches 
up to lOGeV/c. Hence the p+p data were fitted with 
a power law function between 4.5 < pr < 7GeV/c and 
then extrapolated. The systematic uncertainty resulting 
from this extrapolation reaches 20% at lOGeV/c, esti- 
mated from a series of fits, where each time one or more 
randomly selected points are omitted and the remaining 
points are re-fitted. 

So far PHENIX has not measured the p+p spectrum of 
7r° at y/s^ = 39 GeV. Therefore, data from the Fermi- 
lab experiment E706 were used. However, the E706 
acceptance (—1. < |?7| < 0.5) is different from that of 
PHENIX (|?7| < 0.35), and since dN/dr] is not fiat, a pr- 
depcndent correction was applied to the E706 data. This 
correction factor was determined from a pythia simula- 
tion by means of the ratio of yields (normalized per unit 
rapidity) when calculated from the observed yield in the 
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PHENIX and E706 acceptance windows. The system- 
atic uncertainty of the correction is 1-2% at 3 GeV/c but 
reaches 20% at 8GeV/c. 
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FIG. 2: (Color online) Nuclear modification factor (Raa) of 
7r° in Au + Au collisions in most central 0-10% (a) and mid- 
peripheral 40-60% (b). Error bars are the quadratic sum 
of statistical and pT-correlated systematic uncertainties (in- 
cluding systematic uncertainties from the p-|-p-collision refer- 
ence). Boxes around 1 are the quadratic sum of the C-type 
uncertainties combined with the A'coii uncertainties. These 
are fully correlated between different energies. Also shown 
for central collisions are pQCD calculations [31 with regular 
Cronin-effect (solid lines) and with the Cronin-effect reduced 
by a factor of two for all three energies (bands). 

Figure [5] shows the nuclear modification factor of 7r°s 
measured in Au -|- Au collisions at = 39, 62.4 and 

200 GeV (data from Q) as a function of pr for most 
central collisions (a) and 40-60% centrality (b). In the 
most central collisions (0-10%) there is a significant sup- 
pression for all three energies, while in mid-peripheral 
collisions (40-60%) at =39 GeV, Raa is consistent 

with unity above pr > 3GeV/c. 

For 0-10% pQCD calculations [ig, ll7| are also 
shown. The solid curves are calculated with the same 
parametrization that was successful for 200 GeV Au -I- Au 
data (and also 200 GeV Gu + Gu [111). Neither the 62.4, 



nor the 39 GeV data are consistent with the predictions. 
The only qualitative agreement is that the turnover point 
of the Raa curves moves to higher pT with lower collision 
energy as observed in the data. The bands are calculated 
within the same framework but with the Cronin-effect re- 
duced and the energy loss varied by ±10%. The 200 GeV 
data are still well described, the 62.4 GeV data are con- 
sistent within uncertainties, but the 39 GeV Raa, partic- 
ularly the shape, is inconsistent with the corresponding 
band. 

Coupled with the observations that the slopes at high 
Pt become much steeper, but the bulk properties (like el- 
liptic flow, energy density, apparent temperature) change 
only slowly in the collision energy range in question, it is 
quite concievable that hard scattering as a source of par- 



ticles at a given px becomes completely dominant only at 
higher transverse momentum, i.e. jet quenching will be 
"masked" up to higher pt- Note that while the shapes at 
lower Pt are different, at pr >~ 7GeV/c Raa is essen- 
tially the same for the 62.4 and 200 GeV data, irrespec- 
tive of centrality (see also Fig. ^ . While in the 39 GeV 
data Raa also shows a decreasing trend at higher pT, 
unfortunately the px reach of the current data sample 
precludes any conclusion as to what would happen to 
their Raa at even higher px- 
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FIG. 3: (Color online) Nuclear modification factor averaged 
for Pt > 6GeV/c. Uncertainties are shown as error bars 
(statictical) , boxes (sum of pT-uncorrelated and A'coii), boxes 
around one (Type B and C and uncertainties from the p+p- 
coUision reference). 



Figrue [3| shows pr-averaged Raa as a function of the 
number of participants. The averaging was done above 
Pt >6GeV/c. Our first observation is that Raa de- 
creases with increasing centrality even for the lowest- 
energy system. Similarly, as already discussed in the 
context of Fig. O at high enough pT the suppression is 
the same at 62.4 and 200 GeV, at all centralities. This is 
remarkable because the power n of the fit to the spec- 
tra changes approximately by two units from 200 to 
62.4 GeV, so the average momentum loss of the partons 
also has to be different in order to compensate the effect 
of the changing slope. The average momentum loss is 
usually defined by the fractional momentum shift 6pT /pt 
between the corresponding Au -I- Au and TAA-scaled p+p 
spectra as follows. Since the power law tails of the p+p 
and Au -I- Au spectra are similar, they can be fitted si- 
multaneously with the same function and same power n 



f{pT) 



A 



(pt(1 + 5pT/pT)y- 



(2) 



with 5pT being the horizontal shift between the scaled 
p+p and the Au + Au spectra. In panel (a) of Fig. HJ 
the observed fractional momentum shifts are shown for 
central collisions, as a function of the Au -I- Au pT- 

Inclusive single-particle spectra at sufficiently high pT 
and collision energy were predicted to exhibit scaling 
with the variable xt = '2.pT/\fs such that the produc- 
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FIG. 4: (Color online) (a) Fractional momentum shift Spr/pr 
between Au + Au and TAA-scaled p+p data as a function of 
the Au + Au pt- (b) Power rieff of a;T-scaling for p+p and 
Au + Au (minimum bias) at various collision energies. 



tion cross section can be written in a form 1^ 3| 
d^a 1 



E- 



t,(xt,\/s) 



G{xt) 



(3) 



where G{xt) is a universal function and n{xT, \/s) char- 
acterizes the specific process . The scaling power 
ncsixr) between any pair of y/s^ energies is then cal- 
culated as 

log{Yield{xT,^/s^)/Yield{xT,^/s^)) , 
n^s{xT) = , I ^1 ^^ (4) 

In panel (b) of Fig. |4l ncff(a^T) is shown when compar- 
ing invariant- 7r° yields in and Au -I- Au collisions 
at different energies. Both the shape and the magni- 
tude of UefiixT) is similar for the 62.4/200 GeV p-fp 
and Au Au as well as for the 39/200 GeV p-l-p data. 
The rise of n^sixT) at lower xt can be attributed to 
the dominance of soft processes [i^, while at higher xt 
they deviate strongly from leading-twist scaling predic- 
tions [l^, 21 1 . However, the shape of ncff(a;T) in the 
39 and 200 GeV Au + Au comparison is very different 
from all others. It may not even reach its maximum in 
the measured xt range, and its constant rise is similar 
to the rise observed in the Iow-xt (soft) region of the 
other data shown. One possible explanation could be 
that while present, hard scattering is still not the over- 
whelming source of high-p^ tt^s in the currently available 
Pt range in 39 GeV Au -I- Au collisions. 

In summary, the tt*' "pt spectra were measured in 
Au + Au collisions at two different energies, y's„„ — 
39 and 62.4 GeV, and compared to the earlier result 
for ^Js^ — 200 GeV. In all cases the high pT part of 
the invariant yields can be well described with a sin- 
gle power law function. The powers decrease consider- 
ably at lower ^s^-^ , and since the soft processes change 
only slowly with collision energy, jet quenching might be 
"masked" up to higher transverse momenta. The high- 
Vt t!"" yields in Au -I- Au at 62.4 GeV are suppressed. 



and above pr > d GeV/c the data points are compa- 
rable with the 200 GeV results at all centralities. The 
7r° yields in Au + Au at 39 GeV are suppressed in the 
most central collisions, but no suppression is apparent 
in more peripheral collisions. At lower energies, a de- 
creasing momentum shift compensates for the steeper 
slopes at high px, making the -Raa's comparable, in fact, 
identical in the case of 62.4 and 200 GeV. When related 
to 200 GeV, ncff{xT) is similar for 62.4 and 39 GeV p+p 
and 62.4 GeV Au -I- Au, but very different for the 39 GeV 
Au -|- Au data. The new data provided in a wide energy 
range of Au -I- Au collisions will help to constrain energy- 
loss models. 
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